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STELLAR GRANULATION AND INTERFEROMETRY 

A. Chiavassa^ and L. Bigot ^ 


Abstract. Stars are not smooth. Their photosphere is covered by a 
granulation pattern associated with the heat transport by convection. 
The convection-related surface structures have different size, depth, 
and temporal variations with respect to the stellar type. The related 
activity (in addition to other phenomena such as magnetic spots, ro¬ 
tation, dust, etc.) potentially causes bias in stellar parameters de¬ 
termination, radial velocity, chemical abundances determinations, and 
exoplanet transit detections. 

The role of long-baseline interferometric observations in this astro- 
physical context is crucial to characterize the stellar surface dynamics 
and correct the potential biases. In this Chapter, we present how the 
granulation pattern is expected for different kind of stellar types rang¬ 
ing from main sequence to extremely evolved stars of different masses 
and how interferometric techniques help to study their photospheric 
dynamics. 


1 Introduction 

Stellar granulation was observed for the first time on the Sun by Herschel (|1801|) . 
but it was Dawes (118641) who coined the term granules. Eventually, today modern 
telescopes provide direct observations (e.g., Carlsson et al.. l2004|l . The granulation 
pattern is associated with heat transport by convection, on horizontal scales of 
the order of a thousand kilometers (Nordlund et ai, 120091) . On the other hand, 
convection is driven primarily by radiative cooling from a thin thermal boundary 
layer, the layer from which most photons can escape to space. The most prominent 
intensity variations on the solar surface, aside from sunspots and faculae, are 
granules i.e. the bright (hot) areas surrounded by dark (cooler) lanes that tile 
the stellar surface. The horizontal scale on which radiative cooling drives the 
convective motions is linked with the granulation diameter. The bright granules 
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are the locations of upflowing hot plasma, while the dark intergranular lanes are 
the locations of downflowing cool plasma (Nordlund et ai, I1990p . 

Stellar granulation manifests either on resolved observables (e.g., Sun images) 
or unresolved photospheric spectral line in terms of widths, shapes, and strengths. 
In particular, the best observational evidence comes from unresolved spectral lines 
because they combine, in their shapes, important properties such as velocity ampli¬ 
tudes (heavily affecting the line width) and velocity-intensity correlations (Nord¬ 
lund et at. 120091) . Similarly, correlation of velocity and temperature cause char¬ 
acteristic asymmetries of spectral lines as well as net blueshifts (Dravins, 119871 
Gray,[2005l). 

Stellar granulation potentially cause bias in stellar parameters, radial velocity, 
chemical abundances determinations, and exoplanet transit detections. For this 
purpose, large efforts have been made in recent decades to use theoretical mod¬ 
elling of stellar atmospheres to solve multidimensional radiative hydrodynamic 
equations in which convection emerges naturally. These simulations take surface 
inhomogeneities into account (e.g., granulation pattern) and velocity fields and 
are used to predict reliable observables. They cover a substantial portion of the 
Hertzsprung-Russell diagram (Magic et al., 120131 Ludwig et al., I2009|) . includ¬ 
ing the evolutionary phases from the main sequence over the turnoff up to the 
red-giant branch for low-mass stars. These simulations of stellar atmospheres are 
essential for much of contemporary astronomy. 

2 Stellar atmosphere, a pathway to multi-dimensional simulations 

The primary source of information for stellar objects is the light they emit, which 
carries information about the physical conditions at its origin. However, in order 
to interpret the information correctly, one first needs either theoretical or semi- 
empirical models of the atmospheric layers at the surface of stars from where 
the stellar radiation escapes. Many of the observable phenomena occurring on the 
surface of the stars are intimately linked to convection. Moreover, the atmospheric 
temperature stratification in the optically thin region, where the emerging flux 
form, is also affected by the interaction between radiative and convective energy 
transport. To account for all these aspects, it is important to use realistic 3D 
radiative hydrodynamical (RHD) simulations of stellar convection (Nordlund et 
al. 119981 Freytag et al. 120121 Ludwig et al. 120041 Vogler et ai. 120041) . 

RHD simulations numerically solves the time-dependent equations for conser¬ 
vation of mass, momentum, and energy coupled to a realistic treatment of the 
radiative transfer. The simulation domains are of two kinds: (i) box-in-a-star 
simulations (Fig. [U right panel; computational time ranging from few days to few 
weeks depending on the stellar type) cover only a small section of the surface layers 
of the deep convection zone (typically ten pressure scale heights vertically), and 
the numerical box includes about ^ 10 convective cells, which are large enough so 
that the cells are not constrained by the horizontal (cyclic) boundaries, (ii) star- 
in-a-box simulations (Fig. [U left panel; computational time of few months) cover 
the whole convective envelope of the star and have been used to model evolved 
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cool stars like Red Supergiant stars, RSGs (Freytag et al. 120121 Chiavassa et al., 
I2011b|) and Asymptotic Giant Branch stars, AGBs (Freytag et al., I2008II so far. 
Once the surface gravity is lower than logg ~ 1, the box-in-a-star simulations 
become inadequate because of the influence of sphericity becomes important; the 
star-in-a-box global simulations are then needed, but those are highly computer¬ 
time demanding and difficult to run, which is the reason why there are only very 
few simulations available so far. 

RHD simulations employ realistic input physics: (i) updated and reliable equa¬ 
tion of state, which characterize the thermodynamic state of the matter; (ii) the 
radiative transfer includes the most recent continuous absorption coefficients as 
well as atomic and molecular line opacities. The abundances employed in the 
computation of the simulation are the latest chemical composition by Ashland et 
al. ([20091). 

In particular, the radiation transport used in RHD simulations is a crucial 
ingredient either for the resulting atmospheric stratification and the time needed 
for the computation. It cannot be computed monochromatically for all the wave¬ 
length across the spectrum because of the very large computational time needed. 
Therefore The frequency dependance of the radiation field can be calculated with 
two approaches: 

• the gray approximation, which completely ignores the frequency dependence, 
is justified only in the stellar interior and it is inaccurate in the optically thin 
layers. This method is often used for star-in-a-box simulations. 

• the more elaborate scheme accounting for non-gray effects, which is based 
on the idea of opacity binning INordlund. II9821 Nordlund et aZ. ll990|l . The 
basic approximation is the so called multi-group scheme (Ludwig et al.. lI994| 
Vogler et al. I2004[) . In this scheme, the frequencies that reach monochro¬ 
matic optical depth unity within a certain depth range of the model atmo¬ 
sphere will be put into one frequency group. Actual RHD simulations use up 
to 5 and 12 bins {star-in-a-box and box-in-a-star simulations, respectively) 
to reproduce the whole frequency dependence across the full spectrum 

3 Detailed radiative transfer to extract observables 

The extraction of reliable observables require detailed monochromatic calcula¬ 
tions. To overcome the small number of bins in RHD simulations, Chiavassa et 
al. (120091) developed a 3D pure-LTE radiative transfer code called Optim3D that 
make computation between 1000 to 200000 A, with any spectral resolution. The 
code takes into account the Doppler shifts caused by convective motions. The ra¬ 
diative transfer equation is solved monochromatically using extinction coefficients 
pre-tabulated as a function of temperature, density, and wavelength (Fig. [TJ left 
panel). The lookup tables are computed for any kind of chemical compositions 
using the same extensive atomic and molecular opacity data as the latest gen¬ 
eration of MARCS models (Gustafsson et al., I2008|l . Optim3D assumes a zero 
micro-turbulence (i.e., microscale, smaller than the free photon path, non-thermal 






















4 What can the highest angular resolution bring to stellar astrophysics? 



Fig. 1. Left panel: Scheme of the interaction between OptimSD and RHD simulations. 
The RHD simulation (in the star-in-a-box configuration, green to blue star in the middle 
of the numerical cube) thermodynamical variables (temperature, Tijk’, density, pink', and 
velocity, Vink) are used in OptimSD to compute the monochromatic emerging intensity 
along different casting rays. Right panel: RHD simulation in the box-in-a-star 
configuration. 


component of the gas velocity in the region of spectral line formation.) because 
the velocity fields inherent in RHD simulations are expected to self-consistently 
and adequately account for non-thermal Doppler broadening of spectral lines. 

The synergy between RHD simulations and OptimSD allows to tackle several 
astrophysical problems and is based on the systematic production of observables 
like: high and low resolution spectra to study the chemical abundances, velocity 
field, and pulsations of stars as well as their stellar parameters; interferometric 
observables to study the stellar dynamics, stellar parameters, planet transit; and 
images to study the impact of convention on astrometric measurements. 

4 Higher spatial frequencies for the stellar granulation signal 

The size of the convective cells is linked to the pressure scale height at optical- 
depth unity (Freytag et ai, I2001|l . The pressure scale height is defined as 


Hp 


ksTeS 

mg 


(4.1) 


where g is the surface gravity, fcs is the Boltzmann constant and m is the mean 
molecular mass (m = 1.31 x mn = 1-31 x 1.67 x 10“^^ grams, for temperatures 
lower than 10000 K). In the above expression. Tip has the dimension of length. 
Stars with low surface gravity have more diluted atmosphere and lower surface 
temperature while more compact objects are hotter. Fig. [2] show the synthetic 
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Fig. 2. Synthetic images of different stars computed with OptimSD at different wave¬ 
length and from RHD simulations. Top left panel: Procyon with Tea = 6591 ± 43 K, 
logg = 4.01 ± 0.03, radius R = 2.023 ± O.O26f?0, and mass (Chiavassa et al, I2012II . 
The red circle and arrow indicate a particular zone used for the enlargement; Top 
right panel: enlargement of the surface of Procyon. Bottom left panel: AGE star with 
Tea = 2542 ± 220 K, logs = “0.83 ± 0.10, R = 443.5 ± 58.37?©, and M = IMq (Freytag 
et g/.. 120081 Chiavassa et g/. . 12010^11 . Bottom right panel: RSG star with Teff = 3710 ±20 
K, logs — 0.047 ± 0.001, R — 376.7 ± 0.57?©, and M — 6M© (Chiavassa et al.. I2011bll . 


images of stars with different stellar parameters and thus different granulation 
pattern. 

Samadi et al. (|2013L Svensson & Ludwig (I2005P and Chiavassa et al. (|2011ap 
reported a direct correlation between the size of the convective cells and the stel¬ 
lar parameters: larger granules have larger effects on the emerging stellar flux. 
Following this principle, the role of long-baseline interferometric observations is to 
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investigate the dynamics of granulation as a function of stellar parameters: thanks 
to the higher angular resolution, interferometry is the ideal tool for exploring stel¬ 
lar convection in term of cell’s size, intensity contrast and temporal variations. 

4.1 Evolved stars 

Evolved stars, such as RSGs and AGBs, are prime targets for interferometry, 
because of their large diameter, proximity, and high infrared luminosity. These 
stars have effective temperatures lower than ^ 4000 K and surface gravity lower 
than log 5 = 1.0. Their atmospheres, resulting from the RHD simulations, is un¬ 
doubtedly irregular, permeated with structures and dynamics, but depend on the 
wavelength probed (Fig. [2] and Fig. [3]). Moreover, the surface inhomogeneities 
and their temporal evolution induce strong fluctuations on the intensity profiles. 
Fig. 0] displays the comparison between a three-dimensional image representation 
of the intensity in a snapshot of a RSG simulation in the H band with an ad-hoc 
limb-darkened image representation for a RSG star and an uniform disk. The dif¬ 
ference is striking and the resulting RHD surface pattern full of numerous intensity 
spikes. These, though related to the underlying granulation pattern (character¬ 
ized by very large convective cells), are also connected to dynamical effects. The 
emerging intensity depends on (i) the opacity run through the atmosphere (and 
for Teff <4000 K, molecules produce strong absorption in particular in the visible) 
and on (ii) the shocks and waves that dominate at optical depths smaller than~l. 

To search for the effect of the granulation pattern on interferometric visibility 
curves and phases, we use the intensity maps computed with OptimSD (e.g.. 
Fig. [2) and calculated a discrete Fourier transform (FT). We introduced also a 
theoretical spatial frequency scale expressed in units of inverse solar radii (Rq^). 
The conversion between spatial frequencies expressed in the latter scale and in the 
more usual scale of 1/" is given by: 


[11 
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(4.2) 


where v is the spatial frequency, 214.9 is the astronomical unit expressed in solar 
radii, and d is the distance of the observed star. Also useful is the following relation 


B 


A • 0.206265 


(4.3) 


where v is the spatial frequency in arcsec”^ at the observed wavelength A in ^m 
for the baseline B of an interferometer in meters. 

Fig. [5] (bottom panels) shows the visibility curves computed for 36 different 
angles from the intensity map of top panel. Since we do not know the exact posi¬ 
tion of the simulated star with respect to the observations, we applied a statistical 
approach and rotated of 5° the intensity map before computing the FT. Small 
dispersion of the visibility curves (grey thin line) around the uniform and limb- 
darkened disks are noticeable: the measure of the stellar diameter with uniform 
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Fig. 3. Appearance of the stellar surface of an AGB simulation (image from Chiavassa 
and Freytag, I2014bll with respect to the wavelength. The present and future interfero¬ 
metric instrument probing these wavelengths are indicated in lower right corner of every 
image. 


and limb-darkened disk models should take into account this dispersion as an ad¬ 
ditional uncertainty which may be up to 5% at the first null point: for example, 
the uniform disk spans value between 794 and 845 R© for the snapshot in Fig. [5] 
(bottom left panel) in order to fit the minimum and maximum amplitude of visi¬ 
bility fluctuations (Chiavassa et ai. I2009|l . 

Moreover, the visibility dispersion increases clearly with spatial frequency (bottom 
right panel), and visibilities strongly deviate from the parametric cases. It is possi¬ 
ble to characterize the typical size distribution of convective cells on RSGs using in¬ 
terferometric observables. Chiavassa et al. (I2010bp used alpha Ori (RSG star) ob¬ 
servations in H band to tackle this problem and using the following method: after 
















8 What can the highest angular resolution bring to stellar astrophysics? 




Fig. 4. Top panel: three-dimensional image of the emerging intensity from a RHD simu¬ 
lation of RSG star in the H band (image from Chiavassa et a/., [2009}. Bottom left panel: 
from a limb-darkened disk with the stellar parameter of the RHD simulation. Bottom 
right panel: from a uniform disk. 


the computation of the Fourier transform, FT, we obtained I {u, v) = FT [I {x, y)], 
where I (u, v) is a map like in Fig. |SJ The resulting complex number I (u, v) was 
multiplied by low-pass and high-pass filters to extract the information from differ¬ 
ent spatial frequency ranges (corresponding to the visibility lobes). Finally, an in¬ 
verse Fourier transform, FT, was used to obtain the filtered image: /filtered (a;, y) = 


FT 


I (u, v) ■ filter 


Figure [6] (top left panel) shows the filtered images at spatial frequencies, v, 
corresponding to the first lobe of Fig.j^ Since we filtered the signal at high spatial 
frequencies, the image appears blurry and seems to contain only information about 
the stellar radius. However, the top right panel displays the signal related to all 
the frequencies higher than the first lobe: in this image, we clearly do not detect 
the central convective cell of ss 30 mas size (60% of the stellar radius, the size of 
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Fig. 5. Top panel: intensity maps of a RSG star in the H band. Bottom left panel: 
visibility curves (grey thin line) computed for 36 different angles with a step of 5° a 
logarithm scale is applied. The solid black curve is a uniform disk model, with a radius of 
810 Rq. The dashed black line is a partially limb-darkened disk (/ (fj.) = 1 — 0.5 (1 — fj,)) 
with a radius of 822 R©. The dotted dashed line is a fully limb-derkened (7 {p.) — 
1 — 1(1 — p)) with a radius of 830 R©. The stellar parameters of this snapshot are: L = 
98400 L©, R = 836.5 R©, T = 3534 K and log(g) = -0.34. Images from Chiavassa et al., 
[20091 


alpha Ori is about 45 mas at this wavelength). Thus, the first lobe also carries 
information about the presence of large convective cells. Figure [5] (bottom row) 
shows the second lobe with convection-related structures of ^ 10-15 mas, (30% of 
the stellar radius), and the third and fourth lobes with structures smaller than 10 
mas. We can detect convection-related structures of different size using visibility 
measurements at the appropriate spatial frequencies. 
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First lobe, i/ < 30 arcsec ' w > 30 arcsec 
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Fig. 6. Intensity maps filtered at different spatial frequencies qualitatively corresponding 
to the lobes of the visibility curve shown in Fig. [5] The size of alpha Ori is about 45 mas 
at these wavelengths (H band). Images from Chiavassa et al. Il2010bll . 


One last important interferometric observables is the closure phase. The closure 
phase is defined as the phase of the triple product (or bispectrum) of the complex 
visibilities on three baselines, which form a closed loop joining (at least) three 
stations A, B, and C. If the projection of the baseline AB is (ui,wi), that for BC 
is (u 2 , V 2 ), and thus {ui + U 2 , vi + V 2 ) for AC, the closure phase is: 


cj)c{ui,Vi,U 2 ,V 2 ) = arg(I/(ui,Wi) XV{U 2 ,V 2 ) X V*{ui + U 2 ,Vi + V 2 )). 

this procedure removes the atmospheric contribution, leaving the phase infor¬ 
mation of the object visibility unaltered (Monnier. 120071 and l2003jl . Closure phases 
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Fig. 7 . Scatter plot of closure phases in the H band (from the intensity map of top panel 
of Fig. El) of 500 random baseline triangles with a maximum linear extension of 40 m. 
Closure phases are plotted against the longest baseline of the triangle. The upper x-axis 
corresponds to synthetic observations of a RSG simulation at an apparent diameter of 
43.6 mas (which corresponds to a Ori at a distance of 174.3 pc). The axisymmetric case 
is represented by the grey lines. Image from Chiavassa et ai. 120091 
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have the main advantage of being uncorrupted by telescope-specific phase errors, 
including pointing errors, atmospheric piston, and longitudinal dispersion due to 
air and water vapor (Le Bouquin et al. 120121) . Eventually, the information carried 
by the closure phases is intrinsically correlated with the morphology of the stellar 
surface inhomogeneities but depend on the relative orientation of interferometric 
baselines and on the distribution of the telescope station of the diluted aperture 
of an interferometer. 

Figure [3 shows the scatter plot of the closure phase for a RHD simulation of a 
RSG in the H band. The closure phases deviate from zero or ±7r (typical values for 
centrosymmetric objects) already at ~10 m (0.0008 Rq^, scaling the simulation to 
an apparent diameter of 43.6 mas). At higher baselines, it clearly differs from zero 
or ±7r, values that are indicating of a point symmetric brightness distribution. This 
is a clear signature of surface inhomogenities. The characteristic size distribution 
on the stellar surface can also be derived from the closure phase: the contribution 
of small-scale convection-related surface structures increases with frequency. It 
may be very efficient to constrain the level of asymmetry of RSG atmospheres by 
accumulating statistics on closure phase at short and long baselines, since they 
are easily measured to high precision. A small departure from zero immediately 
inters a departure from symmetry (Ghiavassa et al. I2009|) . 

RHD simulations are necessary for a proper quantitative analysis of interfero¬ 
metric observations of the surface of evolved stars beyond the smooth, symmetrical, 
limb-darkened intensity profiles. Frequencies corresponding to the first lobe (ie., 
sensitive to the radius) are characterized by fluctuations as high as 5%, and radii 
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Fig. 8. Snapshot of a RSG simulation convolved to the PSF of a 42 m telescope approx- 
imatively like the European Extremely Large Telescope (for a star at a distance of 152.4 
pc). Image from Chiavassa et al. I2011a1 


determinations can be affected to that extent. The second, third, fourth lobes, and 
so on carry the signature of limb-darkening, and of smaller scale structures due the 
granulation, and are very different from the simple parametric cases. Also closure 
phases largely differ from 0 and ±7r with an increase of the signal with frequency, 
this is due to the departure from circular symmetry. The characteristic size of the 
granulation pattern of, eg., RSGs is composed of convection-related structures of 
different sizes, including small to medium scale granules (5-15 mas) and a large 
convective cell (w30 mas). 

In the last years, RHD simulations managed to explain several interferometric 
observations of RSGs (Chiavassa et al., I2009L I2010b|l and AGBs (Chiavassa et 
ai. 12010^ from the visible to the infrared region. However, the final consistency 
check will be an image reconstruction to directly compare the granulation size 
and shape and the intensity contrast, provided by the planned second generation 
recombiner of the VLTI (MATISSE, GRAVITY, see other contributions is this 
book) and CHARA optical interferometry arrays. First steps in this direction 
has been carried out in Berger et al. (j2012p and Malbet et al. (j2010p where the 
image reconstruction algorithms have been tested using intensity maps from these 
RHD simulations. The European Extremely Large Telescope (E-ELT, planned 
to be operating in 2025) with a mirror size five times larger than a single VLT 
Unit Telescope will be capable of near IR observations of surface details on RSGs 
(Fig. El). 
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Fig. 9. Evolution of the surface inhomogeneities during the stellar evolution: main 
sequence {left panel), red giant {central panel), and AGB {right panel) (the scale is not 
preserved). 


5 Transition from main sequence stars to evolved stars 

The surface inhomogeneities at the surface of the star, i.e. the granulation aspect, 
depend on the effective temperature and even more important the surface gravity. 
The size of the granule is controlled by the stratification and the exponential drop 
of density at the stellar surface (e.g. Nordlund et al. 119901 and Stein et al. 119981) . 
The effective temperature and metallicity play a less important role (Magic and 
Asplund I2014|l . On the Sun, the average granule size is about 1-2 Mm, therefore 
there is at least 10^ granules on the entire surface. The consequence is that even if 
the granule aspect is very irregular, their tiny size with respect to the stellar radius 
makes the appearance of the surface homogenous (Fig.|9]). During the evolution of 
a solar-type star, the radius increases by more than an order of magnitude which 
lower the surface gravity, e.g. logg ^ 2 for the red giant phase and even lower 
than 0 for the AGB phase. The star inflates but at the same time the granule 
size increases to by a factor of 10^ (Magic and Asplund I2014p . The number of 
granules covering the entire surface for a red giant is much smaller than on the 
main sequence phase, with roughly a few thousand : the stellar surface appears 
more and more irregular as the star evolves. 

Clearly, for the latest stage of the stellar evolution the surface of the star is 
very inhomogenous. The operation of computing many Fourier transform (FT) is 
however costly from the numerical point of view. In the process of deriving the 
angular diameter of the star, one needs to compute many time the visibilities of 
the fringes using for example Monte Carlo procedure to derive reliable diameter 
and its error bars. Most of the objects for which we need precise angular diame- 
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Limb darkened scar 
seen by the observer 


Fig. 10. Left panel: Geometry of the system, the star with axial symmetry is represented. 
The surface brightness depends on the position on the stellar disk, r, because of the kimb 
darkening, but does not depend on as long as the granules are very small compared 
with the radius. Right panel: Cut of the star in the plane perpendicular to the left panel. 
The X angle measures the position on the stellar disk 


ters are main sequence stars or early phase of red giants for which we can combine 
asteroseismology and interferometric determination of the diameter to tidily con¬ 
strain the other stellar parameters (e.g., Creevey this book; Bigot et al. I2011l and 
I2015|l . In these cases we can take advantage of the relatively homogenous surface. 
The idea is to simplify the FT using axial symmetry of the star. The steps are 
described below. 

The usual visibility in the (u, u)-plane can be computed using the theorem van 
Citter-Zernike (e.g.. Born and Wolf ll99^ 

F{u, v) = FT [/(x, y)] = j /(x, (5.1) 

This 2D Fourier transform can be written in polar coordinates (r, 9) where r 
is the position of the stellar disk 9 the azimuthal angle. Therefore x = r cos (f and 
y = r sin ip and the surface element can be transformed as dxdy = r dr dip (r being 
the Jacobian of the transformation). The Fourier transforms writes 

F{u,v) = y(5.2) 

It is convenient to write the coordinate (u, u) also in polar coordinates 
Let us, therefore define the following u = k cos (j) and v = k sin if. The Fourier 
transform then writes 
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Fik,^) 


= / 


2'7Tirk{cos (f) cos 


^+"in<^"in‘^)rdrd6» 




(5-3) 

It is important to note that no assumption is made to derive Eq. 15.31 If the 
intensity is now axisymmetric then I (r, ip) = I (r) and we can separate the radial 
and angular integration. The latter is identical to a zeroth-order Bessel function 
of the first kind (Abramowitch and Stegur n~972l) . 


^krcos(v>-<p)^0^2TrJoikr). (5.4) 

In the latter expression, the 4> coordinate disappeared since for any value of (f) we 
can always makeachange of variable so that / = J g27rzfercos(i^')^^/^ 

The Fourier transform of the stellar disk intensity in the case of axial symmetry 
simplifies to a ID integral 



.e/2 

F{k) = 27r / I{r)Jo{kr)rdr, (5.5) 

Jo 

which is the Hankel transform of the disk intensity I (r). Since k is the spatial 
frequency we have k = B/X, where B is the baseline between the two telescopes 
and A is the wavelength of the measurement. The position, r, on the stellar disk 
measures the distance between the disk center r = 0 and the limb r = d/2, where 9 
is the angular diameter of the star. It is common to express since distance is term 
of the angle x made by the vertical (disk center) and the line-of-sight of an observer 
(see Fig. fTOl) . Therefore the position of the stellar disk is simply r = 9 (I — jj?) 
and drjr = —1/(1 — /i^)/id/i where p, = cosx- 

With these definition we car rewrite the visibilities as 


V (X — p?9B/X)fj,dfi (5.6) 

Jo 

which is the usual expression using for fitting angular diameters from the computed 
limb-darkened intensities. 

5.1 Main sequence and giant stars 

As explained in Section [2 the computational domain of RHD simulations of main 
sequence and K giant stars represents only a small portion of the stellar surface. 
To overcome this limitation, and at the same time account for limb-darkening ef¬ 
fects, Chiavassa et al. pOIOcl I20I2112014a|) computed intensity maps for different 
inclinations (d-angles) with respect to the vertical direction and for representa¬ 
tive series (covering several convection turnover) of simulation snapshots and used 
them to tile a spherical surface. The computed values of the d-angle depend on the 
position (longitude and latitude) of the tile on the sphere and are linearly interpo¬ 
lated among the inclination angles. In addition to this, the statistical tile-to-tile 
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fluctuations (i.e., number of granules, shape, and size) is taken in consideration 
by selecting random snapshots within each simulation’s time-series. As a con¬ 
sequence, the simulation assumption of periodic boundary conditions resulted in 
a tiled spherical surface globally displaying an artifactual periodic granulation 
pattern. However, the resulting artificial signal introduced in the interferometric 
observables is less important that the one caused by the inhomogeneities of the 
stellar surface. The final result is an orthographic projection (i.e., stellar disk 
image as seen from a nearby observer) of the tiled spheres (Fig. [2] and ITT]) . 

The extraction of interfererometric observables is done using the stellar disk 
images and the FT method described in the previous sections. However, only 
closure phases are sensitive enough to the effect of stellar granulation to expect a 
possible detection with today interferometers. Chiavassa et al. (|2014a[) computed 
the closure phases for twelve interferometric instruments covering wavelengths 
ranging from visible to infrared, optimizing the observability to allow a broad 
coverage of spatial frequencies. The aim of that work is to present a survey of the 
granulation pattern of stars with different stellar parameters and to evaluate its 
effect on the detection of planet transit. Fig. [TT] displays inhomogeneous brightness 
distributions in the disk images; moreover, the centre-to-limb variations are more 
pronounced in the visible instruments with respect to the infrared ones. This effect 
is explained by different sensitivity of the source (Planck) function at visible and 
at infrared wavelengths. 

Figure [13] (top panel) displays one example from Chiavassa et al. (|2014al) of 
closure phases deviating from the axisymmetric case for MIRC-6 telescope config¬ 
uration case. The scatter plot of closure phase increase particularly in the visible 
wavelength range, where the dispersion is larger (e.g., VEGA, NPOI, and PAVO 
instruments). Depending on the instruments and spatial frequency spanned, the 
departures from symmetry may be large or not. However, it is apparent that the 
convection-related surface structures have a signature on the closure phases. 

Figure [T^ ('bottom panel) shows indeed a direct application to two real targets: 
Beta Com, GOV star with angular diameter of 1.1 mas (Richichi et al., I2005p : 
and Procyon, F5IV with angular diameter of 5.4 mas (Chiavassa et al.. 12012^ . In 
general, all the instruments (except for MATISSE and NPOI, which do not probe 
frequencies larger than the first lobe) show closure phases departures {ijj) of the 
order of few degrees with largest values of the order of ^ 16°. These values are 
much larger than what are the actual instrument’s incertitudes. More in detail: 

• PAVO and VEGA show departures lower than 0.5° already on the 2nd lobe 
(VEGA instrument); 

• AMBER and GRAVITY with values lower than 0.8°; 

• PIONIER in the H band with values of 4.3° (2nd lobe) 6.4° (3rd lobe), and 
PIONIER in the K band with 2.9° (2nd lobe); 

• MIRC with 13.8°, 15.3°, 16.4°, 13.1° for the 3rd, 4th, 5th, 6th lobe, respec¬ 
tively; 
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Fig. 11. Synthetic stellar disk images of different RHD simulations (columns). The 
images correspond to a representative wavelength for each interferometric instruments 
from the visible (top row) to the far infrared (bottom row). The averaged intensity 
(x 10® ergcm“^ s“^ A“^) is reported in the lower left corner of each image. Image from 
Chiavassa et al. Il2014all . 


• CLIMB J band with 6.5°, 0.2°, 5.1° for the 3rd, 4th, 5th lobe, respectively; 

• CLIMB K band with 12.3°, 6.3° for the 3rd, 4th lobe, respectively. 

The actual instruments and telescopes allow, in principle with very good weather 
conditions, the detection of the granulation. The closure phase signal is already 
more pronounced in the infrared for the 2nd lobe and may be detected with very 
good weather and instrumental conditions but it is certainly easier to detect from 
the 3rd lobe on. MIRC instrument with 6 telescope recombination is the most ap¬ 
propriate instrument as it combines good UV coverage and long baselines probed. 

Another important role of long-baseline interferometric observations is planet 
hunting as a complement to the radial velocity and adaptive optics surveys. High 
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Fig. 12. Enlargement of the synthetic stellar disk images of Fig. [TT] for the VEGA 
instrument mounted at CHARA interferometer. Image from Chiavassa et al. Il2014all . 


angular resolution is an ideal tool for exploring separations in the range 1 to 50 mas 
(Le Bouquin et al.. I2012[l . In addition to this, Chiavassa et al. (|2014al) measured 
the contamination of granulation signal on the planetary transit using closure 
phases. To do this, the authors chose three prototypes of planets representing 
different sizes and compositions because the purpose is not to reproduce the exact 
conditions of the planet-star system already detected but to have a statistical 
approach on the interferometric signature for different stellar parameters hosting 
planets with different sizes. Then they created transiting images (Fig. [HI top) and 
computed the resulting closure phases. Eventually, they determined the differences 
between planet-star system and the star alone (Fig. [HI bottom). For all the 
instruments, the absolute difference scales with the size of the planet considered: 
the smaller planet returns smaller differences. In addition to this, the differences 
are larger in the visible wavelengths where the granulation contrast is higher than 
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Fig. 13. Top panel: Scatter plot of closure phases of 20000 random baseline triangles 
(black dots) as a function of the maximum linear extension corresponding to the con¬ 
figuration chosen for each MIRC instrument and for the RHD simulation of the Sun. 
The colored symbols over-plotted display the closure phases for a particular telescope 
configuration chosen and the vertical dashed red lines give the approximative positions 
of the different lobes. Bottom panel: Departure, tp, of closure phases from zero and/or 
±7r (i.e., no axisymmetric case) for RHD simulations with stellar parameters correspond¬ 
ing to two real stars (horizontal axis) and different interferometric instruments (vertical 
axis). 6 lobes are displayed: black for the 1st lobe, violet for the 2nd, light blue for the 
3rd, green for the 4th, yellow for the 5th, and red for the 6th. Only the lobes spanned 
in the UV-planes are plotted. The symbols correspond to different values, in degrees. 
Images from Chiavassa et al. Il2014ap . 
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PIONIER H bond 



Fig. 14. Top: Synthetic stellar disk images in H band (PIONIER instrnment) together 
with three planet transiting phases (black color) for a star like the Sun. The prototypes of 
planet have parameters corresponding to: Kepler-11 f (left), HD 149026 b (central), and 
CoRoT-14 b (right). Bottom: Absolute closure phase differences (in degrees) between 
the star with a transiting planet (from above) and the star alone (Fig. 1111) . The black 
colour correspond to Kepler 11-f, the red to HD 149026 b, and the blue CoRoT 14-b. 
The star symbols connected with solid lines correspond to the planet phase entering in 
the stellar disk, the circle symbols connected with dotted line to the planet at the centre 
of the stellar disk, and the triangles connected with dashed line to the planet exiting the 
stellar disk. Images from Chiavassa et al. (I2014all . 


in the infrared. 

The signature of the transiting planet on the closure phase is mixed with the 
signal due to the convection-related surface structures. The time-scale of granula¬ 
tion depends on the stellar parameters, and varies from minutes or tens of minutes 
for solar type stars and sub-giants, to hours for more evolved red giant stars. If 
the transit is longer that the granulation time-scale (which is the case for most 
of main sequence stars), it is possible to disentangle its signal from convection by 
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observing at particular wavelengths (either in the infrared or in the visible) and 
measuring the closure phases for the star at difference phases of the planetary 
transit. 

For this purpose, it is very important to have a comprehensive knowledge of 
the host star to detect and characterize the orbiting planet, and RHD simulations 
are very important to reach this aim. 


6 Stellar parameters 


Stellar parameters determined with interferometric techniques often consist into 
a simple uniform disk model-fitting that, although unphysical, has the advantage 
of immediately telling if a star is resolved, after which limb-darkening corrections 
are applied. Chiavassa et al. (I2010cll reported the effect on stellar parameter de¬ 
termination of 3D intensity distribution with respect to ID models. This has been 
carried out using synthetic intensity profiles (Fig. (TSl top) from 3D simulations 
with different stellar parameters and metellicities that have been compared to ID 
models with exactly the same parameters. 

The visibility can be derived using V\ {B, 0) from the intensity profile I (A, /r) 
using the Hankel integral: 

Vx{B,Q) = ^j /(A,^) Jo (6.1) 

where A is the wavelength in meters, B is the baseline in meters, 0 is an arbitrary 
angular diameter in radians (2 mas here), Jq the zeroth order of the Bessel function, 
fi = cos{d) (with 6 the angle between the line of sight and the radial direction), 
and A the normalization factor: 


A = 


I (A, /i) /xd/i 


( 6 . 2 ) 


Jo 

For typical red giant stars having 4600 < T^g < 5100 K, the differences in an¬ 
gular diameters vary from about —3.5% to about 1% in the visible, and are roughly 
in the range between —0.5 and —1.5% in the infrared (Table[T|), the corresponding 
change in effective temperature being AT^g/Teg = 1 — •\/ 03 d/ 0 id- 

While the impact of the corrections in Table [T] is usually not dramatic, they 
are not negligible to properly set the zero point of the effective temperature scale 
derived by mean of this fundamental method. In particular in the visible, the (par¬ 
tially) resolved (very) metal-poor stars, it is important to take those corrections 
into account for a correct derivation of their diameters. A second point of concern 
is the reliability of the existing catalogs of calibrator stars for interferometry that 
are based on red giants: the formal high accuracy (a fraction of 1%) based on ID 
models on diameter determination (Merand et al. I2010p may be impacted by the 
correction reported in Table [TJ 

Interferometry is advantageous in that it provides the ability to directly mea¬ 
sure stellar angular diameters and effective temperatures. However, there is a non- 
negligible and still partly unexplored model dependence (Casagrande et oi. ]20I4p . 
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Fig. 15. Top panels: limb-darkening fits (red dashed line) for the RHD azimuthally 
average intensity profile (solid line) for a solar metallicity K giant star (left) and a metal 
poor one (right). The dotted line is the intensity profile computed with the ID model 
having identical stellar parameters, input data, and chemical compositions as the 3D 
simulation. Parametric full (dash-dotted line) and partial limb-darkening (triple dot- 
dashed line) are also shown. Bottom panels: Visibility curves computed with Hankel 
transform with a fixed angular diameter of 2 mas for the above intensity profiles of 
the 3D simulation (solid line), ID (dotted line), and fully limb darkening (dash-dotted 
line). A logarithm scale is used on y-axis. Synthetic visibilities in these plots are not 
realistic near the nulls but are intended only for model-to-model comparison. Images 
from Chiavassa et al. (12010c II . 


More details about the link between stellar parameters, interferometry and 
asteroseismology are reported in Creevey’s contribution to this book. 

7 Additional notes on chromospheric stellar activity with interferome¬ 
try 

The chromosphere (in the classical, stratified, and highly oversimplified view) is an 
intermediate region in the atmosphere of a star, lying above the photosphere and 
below the corona. Chromospheric activity, which encompasses diverse phenom- 
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Table 1. Ratio between the limb-darkened diameters recovered using ID models (0id) 
or 3D simulations ( 03 d) and the corresponding change in effective temperature ATeff for 
the RHD simulations. Table from Chiavassa et al. (12010c I) . 


A [^m] 

[Fe/H] 

log g 

03d/01D 

ATeff [K] 

0.5^ 

0.0 

2.2 

1.003 

-7 

2.14^ 

0.0 

2.2 

0.996 

9 

0.5 

-1.0 

2.2 

0.991 

21 

2.14 

-1.0 

2.2 

0.996 

9 

0.5 

-2.0 

2.2 

0.982 

46 

2.14 

-2.0 

2.2 

0.990 

25 

0.5 

-3.0 

2.2 

1.011 

-28 

2.14 

-3.0 

2.2 

0.989 

28 

0.5 

-3.0 

1.6 

0.965 

82 

2.14 

-3.0 

1.6 

0.984 

37 


“central wavelength of the corresponding visible filter used in Chiavassa et al. 112010c ll 
^central wavelength of the corresponding FLUOR (K band) filter used in Chiavassa et al. 

mwch 


ena that produce emission in excess of that expected from a radiative equilibrium 
atmosphere, is tightly linked to changes in the stellar magnetic field, whether 
periodic or irregular, and is therefore tied to the structure of the subsurface con¬ 
vection zone, the stars rotation, and the regeneration of the magnetic field via a 
self-sustaining dynamo fHall l2008|l . 

Activity signature may manifest in changes in the radiative transfer of lines 
formed very high in the atmosphere of a star, typically for spectral lines sensitive to 
the increase in temperature of the atmosphere. The hydrogen Balmer lines and the 
doublet H&K of the Ca II are used to estimate the strength of the chromospheric 
activity. Mg II lines in the ultraviolet can also help in this direction. 

The advent of a new instruments (VEGA at CHARA), working in near IR and 
visible, opened the possibility to probe and characterize the chromosphere using 
interferometry: for this purpose Ha and Ca II triplet lines have been observed 
using VEGA to extract the signature of stellar activity of a late-type K giant star 
(Berio et al., I2011|) . The authors managed to determine the physical extents of 
the chromosphere of the star measuring the ratio of the radii of the photosphere 
with respect to the chromosphere using the interferometric measurements in the 
Ha and the Ca H infrared triplet line cores. These interferometric measurements 
are unique and crucial to constrain chromospheric models. 
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